skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Leja, Joel"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Free, publicly-accessible full text available November 7, 2025
  2. Abstract We present a population of 11 of the faintest (>25.5 AB mag) short gamma-ray burst (GRB) host galaxies. We model their sparse available observations using the stellar population inference codeProspector-βand develop a novel implementation to incorporate the galaxy mass–radius relation. Assuming these hosts are randomly drawn from the galaxy population and conditioning this draw on their observed flux and size in a few photometric bands, we determine that these hosts have dwarf galaxy stellar masses of 7.0 log ( M * / M ) 9.1 . This is striking as only 14% of short GRB hosts with previous inferred stellar masses hadM*≲ 109M. We further show these short GRBs have smaller physical and host-normalized offsets than the rest of the population, suggesting that the majority of their neutron star (NS) merger progenitors were retained within their hosts. The presumably shallow potentials of these hosts translate to small escape velocities of ∼5.5–80 km s−1, indicative of either low postsupernova systemic velocities or short inspiral times. While short GRBs with identified dwarf host galaxies now comprise ≈14% of the total Swift-detected population, a number are likely missing in the current population, as larger systemic velocities (observed from the Galactic NS population) would result in highly offset short GRBs and less secure host associations. However, the revelation of a population of short GRBs retained in low-mass host galaxies offers a natural explanation for the observedr-process enrichment via NS mergers in Local Group dwarf galaxies, and has implications for gravitational-wave follow-up strategies. 
    more » « less
  3. Abstract A fundamental question in galaxy and black hole evolution remains how galaxies and their supermassive black holes have evolved together over cosmic time. Specifically, it is still unclear how the position of X-ray active galactic nucleus (AGN) host galaxies with respect to the star-forming main sequence (MS) may change with the X-ray luminosity (LX) of the AGN or the stellar mass (M) of the host galaxy. We use data from the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) to probe this issue. XMM-SERVS is covered by the largest medium-depth X-ray survey (with superb supporting multiwavelength data) and thus contains the largest sample to date for study. To ensure consistency, we locally derive the MS from a large reference galaxy sample. In our analysis, we demonstrate that the turnover of the galaxy MS does not allow reliable conclusions to be drawn for high-mass AGNs, and we establish a robust safe regime where the results do not depend upon the choice of MS definition. Under this framework, our results indicate that less massive AGN host galaxies ( log M 9.5 10.5 M ) generally possess enhanced star formation rates compared to their normal-galaxy counterparts while the more massive AGN host galaxies ( log M 10.5 11.5 M ) lie on or below the star-forming MS. Further, we propose an empirical model for how the placement of an AGN with respect to the MS (SFRnorm) evolves as a function of bothMandLX
    more » « less
  4. Abstract We present individual star formation histories (SFHs) of ∼3000 massive galaxies (log(M*/M) > 10.5) from the Large Early Galaxy Astrophysics Census spectroscopic survey at a lookback time of ∼7 billion yr and quantify the population trends leveraging 20 hr deep-integrated spectra of these ∼1800 star-forming and ∼1200 quiescent galaxies at 0.6 <z< 1.0. Essentially all galaxies at this epoch contain stars of age <3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectrophotometric analysis using parametric and nonparametric Bayesian stellar population synthesis modeling tools—BagpipesandProspector—to constrain the median SFHs of this mass complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified ast50andt90, corresponding to the age of the Universe when galaxies formed 50% and 90% of their total stellar mass, although the two methods disagree at the earliest formation times (e.g.,t10). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local Universe—low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star formation (τSF,late=t90–t50) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectrophotometric studies of the more distant Universe. 
    more » « less
  5. Abstract Identifying the sites of r-process nucleosynthesis, a primary mechanism of heavy element production, is a key goal of astrophysics. The discovery of the brightest gamma-ray burst (GRB) to date, GRB 221009A, presented an opportunity to spectroscopically test the idea that r-process elements are produced following the collapse of rapidly rotating massive stars. Here we present James Webb Space Telescope observations of GRB 221009A obtained +168 and +170 rest-frame days after the gamma-ray trigger, and demonstrate that they are well described by a SN 1998bw-like supernova (SN) and power-law afterglow, with no evidence for a component from r-process emission. The SN, with a nickel mass of approximately 0.09 M, is only slightly fainter than the brightness of SN 1998bw at this phase, which indicates that the SN is not an unusual GRB-SN. This demonstrates that the GRB and SN mechanisms are decoupled and that highly energetic GRBs are not likely to produce significant quantities of r-process material, which leaves open the question of whether explosions of massive stars are key sources of r-process elements. Moreover, the host galaxy of GRB 221009A has a very low metallicity of approximately 0.12 Zand strong H2emission at the explosion site, which is consistent with recent star formation, hinting that environmental factors are responsible for its extreme energetics. 
    more » « less
  6. We present a new rest-frame color–color selection method using synthetic us − gs and gs − is, (ugi)s colors to identify star-forming and quiescent galaxies. Our method is similar to the widely used U − V versus V − J (UVJ) diagram. However, UVJ suffers known systematics. Spectroscopic campaigns have shown that UVJ-selected quiescent samples at z ≳ 3 include ∼10%–30% contamination from galaxies with dust-obscured star formation and strong emission lines. Moreover, at z > 3, UVJ colors are extrapolated because the rest-frame band shifts beyond the coverage of the deepest bandpasses at <5 μm (typically Spitzer/IRAC 4.5 μm or future JWST/NIRCam observations). We demonstrate that (ugi)s offers improvements to UVJ at z > 3, and can be applied to galaxies in the JWST era. We apply (ugi)s selection to galaxies at 0.5 < z < 6 from the (observed) 3D-HST and UltraVISTA catalogs, and to the (simulated) JAGUAR catalogs. We show that extrapolation can affect (V − J)0 color by up to 1 mag, but changes $${({g}_{s}-{i}_{s})}_{0}$$ color by ≤0.2 mag, even at z ≃ 6. While (ugi)s-selected quiescent samples are comparable to UVJ in completeness (both achieve ∼85%–90% at z = 3–3.5), (ugi)s reduces contamination in quiescent samples by nearly a factor of 2, from ≃35% to ≃17% at z = 3, and from ≃60% to ≃33% at z = 6. This leads to improvements in the true-to-false-positive ratio (TP/FP), where we find TP/FP ≳2.2 for (ugi)s at z ≃ 3.5 − 6, compared to TP/FP < 1 for UVJ-selected samples. This indicates that contaminants will outnumber true quiescent galaxies in UVJ at these redshifts, while (ugi)s will provide higher-fidelity samples. 
    more » « less
  7. Abstract One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the simba cosmological simulation with the prospector SED fitting code. We consider the impact of the commonly assumed uniform screen model and introduce a new nonuniform screen model parameterized by the fraction of unobscured stellar light. This nonuniform screen model allows for a nonzero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a nonuniform screen model, with the median offset in UV attenuation decreasing from −0.30 dex with a uniform screen model to −0.17 dex with the nonuniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the nonuniform screen model, decreasing the SFR offset on average by 0.12 dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations. 
    more » « less
  8. Abstract Poststarburst galaxies (PSBs) are young quiescent galaxies that have recently experienced a rapid decrease in star formation, allowing us to probe the fast-quenching period of galaxy evolution. In this work, we obtained Hubble Space Telescope (HST)/WFC3 F110W imaging to measure the sizes of 171 massive ( log ( M * / M ) 11 ) spectroscopically identified PSBs at 1 <z1.3 selected from the DESI Survey Validation luminous red galaxy sample. This statistical sample constitutes an order of magnitude increase from the ∼20 PSBs with space-based imaging and deep spectroscopy. We perform structural fitting of the target galaxies withpysersicand compare them to quiescent and star-forming galaxies in the 3D-HST survey. We find that these PSBs are more compact than the general population of quiescent galaxies, lying systematically ∼0.1 dex below the established size–mass relation. However, their central surface mass densities are similar to those of their quiescent counterparts ( log ( Σ 1 kpc / ( M kpc 2 ) ) 10.1 ). These findings are easily reconciled by later ex situ growth via minor mergers or a slight progenitor bias. These PSBs are round in projection (b/amedian∼ 0.8), suggesting that they are primarily spheroids, not disks, in 3D. We find no correlation between the time since quenching and light-weighted PSB sizes or central densities. This disfavors apparent structural growth due to the fading of centralized starbursts in this galaxy population. Instead, we posit that the fast quenching of massive galaxies at this epoch occurs preferentially in galaxies with preexisting compact structures. 
    more » « less
  9. Abstract W-CDF-S, ELAIS-S1, and XMM-LSS will be three Deep-Drilling Fields (DDFs) of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), but their extensive multiwavelength data have not been fully utilized as done in the COSMOS field, another LSST DDF. To prepare for future science, we fit source spectral energy distributions (SEDs) from X-ray to far-infrared in these three fields mainly to derive galaxy stellar masses and star formation rates. We use CIGALE v2022.0, a code that has been regularly developed and evaluated, for the SED fitting. Our catalog includes 0.8 million sources covering 4.9 deg 2 in W-CDF-S, 0.8 million sources covering 3.4 deg 2 in ELAIS-S1, and 1.2 million sources covering 4.9 deg 2 in XMM-LSS. Besides fitting normal galaxies, we also select candidates that may host active galactic nuclei (AGNs) or are experiencing recent star formation variations and use models specifically designed for these sources to fit their SEDs; this increases the utility of our catalog for various projects in the future. We calibrate our measurements by comparison with those in well-studied smaller regions and briefly discuss the implications of our results. We also perform detailed tests of the completeness and purity of SED-selected AGNs. Our data can be retrieved from a public website. 
    more » « less